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Abstract 

In this paper we give an extension of the Birkhoff-Lewis theorem to some 
semihnear PDEs. Accordingly we prove existence of infinitely many periodic orbits 
with large period accumulating at the origin. Such periodic orbits bifurcate from 
resonant finite dimensional invariant tori of the fourth order normal form of the 
system. Besides standard nonresonance and nondegeneracy assumptions, our main 
result is obtained assuming a regularizing property of the nonlinearity. We apply 
our main theorem to a semilinear beam equation and to a nonlinear Schrodinger 
equation with smoothing nonlinearity. 

1 Introduction 

In 1934 Birkhoff and Lewis [BL34] (see also [Lew34, Mos77]) proved their celebrated 
theorem on existence of periodic orbits with large period close to elliptic equilibria of 
Hamiltonian systems^. Here we give a generalization of their result to some semilinear 
Hamiltonian PDEs. 

Birkhoff-Lewis procedure consists in putting the system in fourth order (Birkhoff) 
normal form, namely in the form 

H = Ho + G, + R,, Ho:=Y.^j^^^^, (1) 

i=i 

where G4 is a homogeneous polynomial of degree 2 in the actions Ij := (p^ + g|)/2 and 
is a remainder having a zero of fifth order at the origin. Then system (1) is a perturbation 
of the integrable system Hq + G4. Under a nondegeneracy condition (that also plays a 

*Dipartimento di Matematica, Via Saldini 50, 20133, Milano, Italy. Supported by GNFM. 

tSISSA, Via Beirut 2-4, 34014, Trieste, Italy. Supported by MURST, Variational Methods and Non- 
linear Differential Equations. 

^Actually [BL34] considers the neighborhood of an elliptic, non-constant, periodic orbit, but the 
scheme is essentially the same for elliptic equilibria. 
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fundamental role in KAM theory) the action to frequency map of this integrable system is 
one to one, and therefore there exist infinitely many resonant tori on which the motion is 
periodic. The question is: Do some of these periodic orbits persist under the perturbation 
due to the term R^l Birkhoff-Lewis used the implicit function theorem and a topological 
argument to prove that there exists a sequence of resonant tori accumulating at the origin 
with the property that at least two periodic orbits bifurcate from each one of them. 

In order to extend this result to infinite dimensional systems describing Hamiltonian 
PDEs one meets two difficulties: the first one is the generalization of Birkhoff normal 
form to PDEs and the second one is the appearing of a small denominator problem. 

Here wc decide to work in a way which is as straightforward as possible, so, instead of 
considering the standard Birkhoff normal form of the system, whose extension to PDEs 
is not completely understood at present^, wc consider its "seminormal form", namely 
the kind of normal form employed to construct lower dimensional tori. Precisely, having 
fixed a positive n, we split the phase variables in two groups, namely the variables with 
index smaller than n and the variables with index larger than n. We will denote by z the 
whole set of variables with index larger than n. We construct a canonical transformation 
putting the system in the form 

Ho + G + G + K (2) 

where G depends only on the actions, G is at least cubic in the variables z with index 
larger than n, and K has a zero of sixth order at the origin. The interest of such a semi- 
normal form is that the normalized system Hq + G + G has the invariant 2n-dimensional 
manifold z = which is filled by n-dimensional invariant tori. 

Under a nondegeneracy condition, the frequencies of the flow in such tori cover an 
open subset of R". We concentrate on the resonant tori filled by periodic orbits and 
we prove that at least n geometrically distinct periodic orbits of each torus survive the 
perturbation due to the term K. Since the orbits bifurcate from lower dimensional tori we 
have to impose a further nondegeneracy condition in order to avoid resonances between 
the frequency of the periodic orbit and the frequencies of the transversal oscillations. 

The proof is based on a variational Lyapunov Schmidt reduction similar to that 
employed in [BBV03] and inspired to [ACE87]. It turns out that in the present case 
the range equation involves small denominators. To solve the corresponding problem we 
use an approach similar to that of [BamOO]. In particular wc impose a strong condition 
on the small denominators and we show that, if the vector field of the nonlinear ity is 
smoothing, then the range equation can be solved by the contraction mapping principle. 
Next, the kernel equation is solved by noting that it is the Euler-Lagrange equation of 
the action functional restricted to the solutions of the range equation. The restricted 
functional turns out to be defined on T" and so existence and multiplicity of solutions 
(critical points) follows by the classical Lusternik-Schnirelmann theory. 

Finally we apply the general theorem to the nonlinear beam equation 

utt + u^xxx + mu^ f{u) . (3) 
^See however the recent works [Bam03,BG03]. 
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with Dirichlct boundary conditions on a segment. We consider m as a parameter varying 
in the segment [0, L] and we show that the assumptions of the abstract theorem are 
fulfilled provided one excludes from the interval a finite number of values of m. As a 
second application we will deal with a nonlinear Schrodinger equation with a smoothing 
nonlinearity of the type considered in [Pos02]. 

We recall that families of periodic solutions to Hamiltonian PDEs have been con- 
structed by many authors (see e.g. [Kuk93m, Way90, Pos96b, CW93, Bou95, BamOO]). 
The main difference is that the periodic orbits of the above quoted papers are a contin- 
uation of the linear normal modes to the nonlinear system. In particular their period is 
close to one of the periods of the linearized system. Moreover (except in the resonant 
case, see [BP01,BB03]) each periodic solution involve only one of the linear oscillators'^. 

On the contrary, the periodic orbits constructed in the present paper are the shadows 
of resonant tori; they are a purely nonlinear phenomenon, have long period, and moreover 
each periodic motion involves n linear oscillators that oscillate with amplitudes of the 
same order of magnitude. 

Acknowledgements: The authors thanks L. Biasco, P. Bolle and M. Procesi for inter- 
esting discussions. 



2 Main result 

Consider a real Hamiltonian system with real^ Hamiltonian function 

H{z, z) = Y^ LOjZjZj + P{z, z) = Ho + P (4) 

where P has a zero of third order at the origin and the symplectic structure is given by 
iJ2j dzj A dzj. Here z and z are considered as independent variables. Often we will write 
only the equation for z since the equation for z is obtained by complex conjugation. The 
formal Hamiltonian vector field of the system is Xh{z,z) :— (i|^, ^-^d therefore 

the equations of motion have the form 

.dP ^ . _ .dP 
Zj = lujjZj + 1— , Zj = -lujjZj - 1—. (5) 

Define the complex Hilbert space 

H'^'^iC) := {w = {w^,W2, . . .) e C- I \\w\\l, := ^ k.lY'e^^" < oo}. 

We fix s > and a > and we will study the system in the phase space 

Va,s ■.^W''{C)XW{C)3{Z,Z) . 

Fix any finite integer n > 2 and denote uj :— {cui, . . . ,a;„), fl := {LJn+i,uin+2, ■ ■ •)• 
We assume 



^in the sense that all the other ones have a much smaller amplitude of oscillation. 

■^i.e. if z is actually the complex conjugate to z then the Hamiltonian H takes real values. 
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(A) The frequencies grow at least linearly at infinity, namely there exists a > and 
di > 1 such that 

(NR) For any A; e Z", Z e Z°° with \l\ < 2 and \k\ + \l\ < 5, one has 

u-k + n-ly^O. (6) 

(S) There exists a neighbourhood of the origin U <Z Vas and c? > such that Xp e 

Remark 2.1 In applications to PDEs, property (S) is usually a consequence of the 
smoothness of the Nemitsky operator defined by the nonlinear part of the equation. In or- 
der to ensure (S) one has usually to restrict to the case where the functions with Fourier 
coefficients in Va,s is an algebra (with the product of convolution between sequences). 
This imposes some limitations to the choice of the indices a,s. 

Proposition 2.1 Assume (A,NR,S). There exists a real analytic, symplectic change of 
variables T defined in some neighborhood W C Va,s of the origin, transforming the 
Hamiltonian H in seminormal form up to order six, namely into 

HoT = n^Ho + G + G + K (7) 
with ^ 

^ ~ 2 ^-^ ^iil^il l^il ' 

min(i j)<n 

Gij = Gji, G = 0(11^11^ J where z := (z„+i,2;„+2, . . .) and K = 0{\\z\\l^^). Moreover 

X^,X^,XKeC^{U',Va,s+d) , \\z-T{z)\\a,s+d<C\\z\\l, . (8) 

We defer the proof of this proposition to the Appendix. 

The interest of such a seminormal form is that the system obtained by neglecting 
the reminder K has the invariant manifold z = on which the system is integrable. 

As a variant with respect to the standard finite dimensional procedure we have left 
the third order term G but normalized the system up to order six (instead of five) . This 
is needed in lemma 3.2. 

We rewrite the Hamiltonian Ti in the form 

n:=uJ-I + n-Z + l-AI-I+{BI,Z) + G + K (9) 

where / := (l-^ip, ■ ■ ■ , l-^np), ^ '■— (l-^n+ip, |-2^n+2P, ■ ■ ■) are the actions, A is the n x n 
matrix 

A = {Gij)i<ij<n (10) 

and B is the oo x n matrix 

B = (G^ii)l<i<n<i (11) 
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Remark 2.2 Due to (8,9), one has \{BI)j\ < C\I\j~'^ for a suitable C. Indeed, since 
Xq maps Va,s to Va,s+d, the operator Zj i— >• {BI)jZj maps Va,s to Va,s+d (ind therefore its 
eigenvalues {BI)j must fulfill the above property. 



Introduce action angle variables for the first n modes by Zj = \zj\e^'^^ = ■\jlje^'''^ for 
J = 1, ...,n. 

Perform the rescaling Ij —>■ rflj, (pj (f)j for j = 1, ... ,n, zj — > rjzj, Zj — > rjZj for 
j > n + 1 and divide the Hamiltonian by rj'^. We get 

n{I, (P,z,%)=LU-I + n-Z + r]Gr, + r]^ {^AI ■ I + {BI, Z)) + r^^Kr,. (12) 

where = Odl-^lla.J and ir^(z) = 0(||^||^ J. We will still denote by Va,s = R"x T"x 
7i"'*x Ti"'"* the phase space. 

We will find periodic solutions of the Hamiltonian system (12) close to periodic 
solutions of the integrable Hamiltonian system 

7 = 0, 4)^u + rf{AI + B'^Z), zj ^ i(nj + rj'^{B I) j^j, j>n + l (13) 

in which and K^j are neglected. The manifold {z = 0} is invariant for the Hamiltonian 
system (13) and it is completely filled up by the invariant tori 

T(7o) := {/ = Io, <pe z = 0} 

on which the motion is linear with frequencies 

CO = u){Io) := a; + rfAI^. 

Such a torus is linearly stable and the frequencies of small oscillation about the torus 
T(7o) are the "shifted elliptic frequencies" , namely 

^•(7o):=(l] + ^^fi/o),. (14) 

If all the cD's are integer multiples of a single frequency, namely if 

a) := a; + rfAI^ = ^27rA; e ^27rZ", (15) 

then T(7o) is a completely resonant torus, supporting the family of T-periodic motions 

P:={7(t)=7o, 0(t)=0o + c5i, i(t)=0}. (16) 

The whole family V will not persist in the dynamics of the complete Hamiltonian system 
(12). We will show that, under suitable assumptions, at least n geometrically distinct 
T-periodic solutions persist. More precisely, we will show that this happens for 77 small 
enough and for any choice of Iq and T with 

WhW <C, 1 < T < 4 (17) 
where C is independent of rj, fulfilling 
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(HI) Equation (15) holds. 

(H2) There exist S > and r < d such that 



l^jT -2Trl\> —, VZeZ, Vj>n + 1. 



(18) 



Proposition 2.2 Fix r > 1. Assume (A), detA ^ and 

Clj ■= (n - BA-^uj) . 7^ , Vj > n + 1, 



(19) 



then, for any rj > 0, and almost any T fulfilling (17) there exists Iq such that (H1,H2) 
hold. 

Proof. Fix ?7, we define Iq := Io{T) as a function of T so that that (15) is identically 
satisfied. Then we find T so that the non resonance property (18) holds. Fix rj and define, 



lo := /o(T):= 
k := k{T) = 



27r 

ujT 



27r. 



27r 



(20) 
(21) 



where [{xi, . . . ,Xn)] '■= {[xi], ■ ■ ■ , [.t„]) and [x] G Z denotes the integer part of a; G R. 
With the choice (20), (21), ujT + Tr]^AIo = 27rk, and Iq is of order 1 since Tr]^ > 1. 

We come to the nonresonance property (18). To study it remark that the function 
T — > [uT/2t:\ is piecewise constant. Hence, for any Tq G [r]^^ , 21]^'^) there exists an 
interval Iq = (Tq — a,To + 6) C [rj^^ , 2ri^'^] such that [c<jT/27r] := /cq is constant for 
T & Iq. Moreover, the union of such intervals cover the whole set of values in which we 
are interested. We will construct a subset of full measure of Iq, in which condition (H2) 
is fulfilled. 

So, for fixed j, I consider the set 



Bji{r, S) := {T G Jo : I^T - 27r/| < S/f} . 



(22) 



Remark that 



so that, in Iq 



n.T = n.T + (2nBA-^ 



2n 



— {Q,T-2nl)=Q,. 



dT 

By (A) and remark 2.2 there exists C such that 



>Cj 



di 



>Cj 



di 



(23) 
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Then Bji is an interval with length \Bji\ controlled by 

\Bji\ < 2 



(24) 



Fix j and estimate the number of I for which the set Bji is (possibly) non empty. First 
remark that, due to (A), one has that, as T varies in Iq, the quantity CljT varies in a 
segment of length smaller than Cj'^^ , with a suitable C. This means that there are at 
most Cj*^^ values of I which fall in such an interval (with redefined C). So, one has 



I 



< 



(25) 



Thus, provided r > 1 as we assumed, one has that 



< cs. 



(26) 



By this estimate, the intersection over S of such sets has zero measure, which is the 
thesis. □ 



Theorem 2.1 Let T and Iq fulfill (17) and (H1,H2), then, provided rj is small enough, 
there exist n geometrically distinct periodic orbits of the Hamiltonian system 7i cfr. (12) 
with period T which are rf close in Va,s to the torus T{Iq). 



Going back to the original system one has 



Corollary 2.1 Consider the Hamiltonian system (5) and fix a positive n. Assume that 
(A,NR,S) hold, that detA ^ (cf. (10)) and that Clj 7^ for all j > n + 1 (cf. (19)). 
Finally assume d> 1. 

Then, for any positive 77 <^ 1 there exist at least n distinct periodic orbits z^^\t), z^'^\t) 
with the following properties: 

• ||2;W(t)||a,s < Ct) for I = 1, ...,n and t e R; 

• ||n>„2;'^'^(t)||a^s < Crj^ for I = 1, ...,n and t G R; here n>„ is the projector on the 
modes with index larger than n; 

• The period T of z^'"^ does not depend on I and fulfills r]~^ <T< 2r]~^ . 



Remcirk 2.3 If the integer numbers {ki,...,kn) — ujT /2t: are relatively prime then T 

is the minimal period of the periodic solutions z^^\ Indeed z'^^^ are T -periodic functions 
close to the functions defined in (16) which have minimal period T. 
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3 Proof of theorem 2.1 



Since the problem is Hamiltonian, any periodic solution of the system is a critical point 
of the action functional 

S(1, 0, 5, t) = r ll-^ + i V zjij - n(1, 0, z, t) I dt (27) 

V j>n+l J 

in the space of T-periodic, Po^g-valued functions. Here Tt is given by (12). 
We look for a periodic solution ( :— (0, /, z, z) of the form 

(P{t) = (l)o + ojt + 'ilj{t), I{t)^ Io + J{t), (28) 

where {ip, J, z, 'z) are periodic functions of period T taking values in the covering space 
R" X R" X H"''" X K"-'-'' of Va,s (that for simplicity wiU still be denoted by Va,s)- Hence 
{ip, J, z) must satisfy (in the sequel for simplicity of notation we will only consider the 
equation for z) 

V^S{C) = ^ j = RdC) (29) 

ViS{O = ^ ij-r]'AJ = Ri{C) (30) 

V^,5(C) = ^ z,-in,z, = {R^),{C) (31) 

where 

R^{C) := -rj'^d^K^{Io + J,(l)o + (^t + 'ilj,z) -rid^Gr, 

< RiiO := r]^B'^Z + r]^diKr,{Io + J,(f)o + Ojt + jjj,z)+ridiGr, (32) 
, {Rz)jiC) = iV^{BJ)jZj + irjdsjGr, + \r]'^d^.K,^{lQ + J,(f)o + u;t + ip, z). 

Remark that, since one expects J, ip and z to be small (they will turn out to be of 
order 77^), and by proposition 2.1 one has d^Gr,{C) = C>(||i||^ J, d^Gr,{C) = 0(||i||^ J, 
diGjj{() — 0(11^11^ 3) the r.h.s. of (29,30,31) is actually small with respect to the period 

T < 2?7~^, see lemma 3.2. 

Define the Hilbert space Hp ((0, T); Va,s) of the T-periodic Va,s valued periodic func- 
tion of class H^. In order to simplify notations we will denote this space by Hp^. Denote 

\J\h,T ■■= ^ [ \J? dt , l^li.,^ := ^ 1^1^ dt, (33) 
\Mh,T,a,s ■■= ^JjHmldt, (34) 

||C||L2,T,a,s \J\l^T + I'^Il^T + \\w\\L^,T,a,s ■ (35) 

We will endow Hp ((0, T);Va,s) = Hp^^ with the norm 

\L^T,a,s + T\\C\\L^T,a,s. (36) 
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Remctrk 3.1 With this choice one has 



IICWIk.<C||Clka,., WeR 

with a constant independent ofT. Therefore, with this choice of the norm, the space Hp^ 
is a 'Banach algebra', and the T. a. s norm of the product of any component of a vector 
( with any component of a vector (' is bounded by C\\(\\T,a,s\\C'\\T,a,s with a constant C 
independent ofT. 

We will consider the system (29,30,31) as a functional equation in Hpg. 

Remark 3.2 As a consequence of (8) and remark 3.1, the map C, ^ R{C) '■— {R<j>{C)j 
RiiC): RziO) «^ « C'^ rnap from H\,^^ to H\^^^. 

We are going to use the method of Lyapunov-Schmidt decomposition in order to 
solve (29,30,31). To this end remark that the kernel of the linear operator C at l.h.s. 
of (29,30,31) is given by (0,0,0) with constant G T". The range of C is the space 
of the functions C = ijP^ J,z) with ip^t) having zero mean value. So, there is a natural 
decomposition of Hp^ in Range+Kernel. Explicitly we write 

C = (0 + ^, J, 5) = (^, J, z) + (0, 0, 0) = + 

with '0 having zero mean value and being constant. Then we fix 0, take the projection 
of the system (29,30,31) on the range and solve it. The solution is a function C-r(^) 4>)- 
Finally we insert this function in the variational principle in order to find critical points 
of S. 



3.1 The Range equation 

The range equation has the form 

j = i?^(C) - (i?^(C)) 

iP-r)lAJ = RiiC) (37) 

zj-iQjZj = (Rz)j{C), 

where {R^{()) :— (l/T) R^{()dt. We look for its solution in the range, namely in the 
space 

Hp,s C H}>^ 

of the functions Cr= {'4^i J-, ^) with ip having zero average. 

First of all we analyze the linear problem defined by the left hand side of (37). A 
"small denominator problem" appears since inverting this linear system the denominators 
(ijT — 27r/, J > n + 1, / e Z are present. So, define the linear operator 

£(■0, J, z) = £(^R := ( j, ip — ffAJ., Wj — i^jWj) 

and study the equation 

CCR = {i^J,w) (38) 

with (^5, J, w) e -f^Ps+T given. 
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Lemma 3.1 Assume (H2). If 

Cr = (V^, e -ffp,s+r i.e. in i?p,^+^ with / ^{t) dt = 0, 

t/ 

then the equation (38) has a unique solution 

Cr= {'tp,J,w) eHp^^. 
Moreover, for T e {r)~^, 2r)~'^) and a constant C :— C{S) 

C ~ 

||Cr||t,o,s < II C-R II T,a,s+T. 

Proof. Since A is symmetric and invertible it has an orthonormal basis of eigenvectors 
ei, . . . , Cn with eigenvalues Ai, . . . , A„. In these coordinates J(t) = J2k=i Jk(t)^k-i i^if) — 
YIk=i '4'k{t)Gk and the solution (r of (38) with ipo = has Fourier coefficients 

Jkl = -77— r for / 7^ 0, J^o = 5-—, 

, rr^Jkl + V^Jkl^k r l/n 

i^ki = T — for / ^ 0, 

IZTT/ 



and, for j' > n + 1 



Twji 



i(27r/ - njT) 



We find 

l^li^r = E \J\l, = i: + i: (^)' < §i\J\h.T + Cni-lLr- (39) 

ki k 'I k,l^O '/ 

A similar estimate for |V'|i,2,t holds. Moreover 

|V'|l2,t < vVIl^t + li^lL^T < C{\J\l2,t + li^lL^r), (40) 
using (39). Finally, the solution w = {'Wj)j>n+i of (38) is 



^^.(t)=^__i:f^e^(WT)^* 



where Wj{t) = Eiez ^jic'^^'^/^)'*. From (H2) we get 

lklU2,T,a,s < Cj\\w\\L^T,a,s+T, ||^i'||L2 ,T,a,s < C j\\w\\l2 ,T,a,s+r- (41) 

By (39), (40) and (41) the last estimate of the lemma follows. □ 
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Thus L ^ defines a Hnear bounded operator L : i7p^^^ — > iJp^. 

In order to find a solution ^p = (■0, J, z) of the range equation it is sufficient to find 
a fixed point of 

Ci? = <f(Ci?):=^(iV(Cij;0)) (42) 
in the space -ffp^, where N :— N{(ji; 0) denotes the right hand side of (37). 

Lemma 3.2 Assume d > t. Then there exists a constant C sufficiently large such that 
V77 <^ 1 the map ^ is a contraction of a ball of radius Crf'. 

Proof. Consider a C_r G ifp^ with ||Ci?||r,a,s < P with some positive (small) p. Since Hp^ 
is an algebra with constants independent of T (c.f. remark 3.1), one has, by (32), 

\\N{CR)\\T,a,s+d<C{n' + r^p') 

with a suitable C. Therefore, by lemma 3.1 one has 

\mR)\\T,a,s < \mR)\\T,a,s+d-r < + ^) 

which is smaller than p provided Cirf + p^ /rj) < p, which is implied e.g. by p = 207}^ 
and rj small enough. 

Similarly one estimates the Lipschitz constant of $ by the norm of its differential. 
Such a differential is bounded in a ball of radius p by C{rj^ + p/rj), from which the thesis 
follows. □ 

Corollary 3.1 There exists a unique smooth function T" 3 (f) Cr{4>,v) ^ -f^p,s solving 
(37) and fulfilling 

\\CR{<P,V)\\T,a,s<Cr]\ 



3.2 The kernel equation 

The geometric interpretation of the construction of the previous subsection is that we 
have found a submanifold T" = (0, C_r(0, v)) <^ ^Ps^ diffcomorphic to an n dimensional 
torus, on which the partial derivative of the action functional S with respect to the 
variables (r vanishes. Consider now the restriction Sn of S to T". At a critical point of 
Sn all the partial derivatives of the complete functional S vanish, and therefore we can 
conclude that such point is critical also for the non restricted functional. 

By standard Lusternik-Schnirelmann theory there exist at least n geometrically 
distinct T-periodic solutions, i.e. solutions not obtained one from each other simply by 
time-translations. Indeed, restrict Sn to the plane E := [uj]-^ orthogonal to the periodic 
flow uj = {1/T)27ik with A; e Z". The set Z" n is a lattice of E, and hence deflnes 
a functional Sn\r on the quotient space F :— E / {Z"- H E) r-j T"~^ 

Due to the invariance of Sn with respect to the time shift, a critical point of Sn\T is 
also a critical point of S'„ : T** ^ R. By the Lusternik-Schnirelman category theory since 
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catr = catT"^^ = n, we can define the n min-max critical values Ci < C2 < . . . < c„ 
for Sn\r- If the critical levels Cj are all distinct, the corresponding T-periodic solutions 
are geometrically distinct, since their actions q are all different. On the other hand, if 
some min-max critical level q coincide, then, by the Lusternik-Schnirelmann theory, Sn\r 
possesses infinitely many critical points. However not all the corresponding T-periodic 
solutions are necessarily geometrically distinct, since two different critical points could 
belong to the same orbit. Nevertheless, since a periodic solution can cross T at most 
a finite number of times, the existence of infinitely many geometrically distinct orbits 
follows. For further details see [BBV03]. 

This concludes the proof of the theorem 2.1. 



4 Applications 

4.1 The nonlinear beam equation 

Consider the beam equation 

utt + u^xxx + mu = f{u) (43) 
subject to hinged boundary conditions 

m(0, t) = u{7r, t) = u^^{0, t) = u^^{n, t) = 0, (44) 
where the nonlinearity f{u) is a real analytic, odd function of the form 

f{u) ^au^ + J2 fku'', a 7^ 0. 

fe>5 

The beam equation (43) is a Hamiltonian PDE with associated Hamiltonian 



where g{u) := Jq f{s) ds is a primitive of /. 
Write the system in first order form 

\v = -Ua:xxx -mu + f{u) 

The standard phase space"'' for (45) is J^s '■— Hq x 3 {u,v), where Hq is the space 

of the functions which extend to skew symmetric H"^ periodic functions over [— tt.tt]. 
Note that = {u{x) = Z)j>i sin(ja;) | J2j>i I'^^jl^J^^ < +00}. It is then immediate 
to reahze that, due to the regularity and the skewsymmetry of the vector field of the 



^An equivalent definition makes use of the so called compatibility conditions required for the smooth- 
ness of solutions of second order equations with Dirichlet boundary conditions, see e.g. [Bre93], theorem 
X.8. 
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nonlinear part, / defines a smoothing operator, namely a smooth map from J^g to 
provided s > 1. 

Here we are also interested in spaces of analytic functions, namely functions whose 
Fourier coefficients belong to with some positive a. It is easy to see that the smooth- 
ing property of the nonlinearity holds also for these spaces. 

Introduce coordinates q— {qi,q2, ■■■,), p — (pi,p2, ■ ■ ■ , ) through the relations 



where = ^J2/Trsm{jx) and 

Remark also that 

Passing to complex coordinates 



u;^ — + m. 



qj + ipj ^ qj-ipj 



(46) 



v/2 ' V2 ' 

the Hamiltonian takes the form (4) and the nonlinearity fulfills (S) with s > 1 a suitable 
a, depending on the anayticity strip of /, and d = 2 (for more details see [Pos96b]- 
[GY03]). 

In order to verify the nonresonance property we use m as a parameter belonging to 
the set [0, L] with an arbitrary L. 

Lemma 4.1 There exists a finite set A C [0, L] such that, ifm e [0, L\\A then condition 
(NR) holds. 

Proof. First remark that, due to the growth property of the frequencies there is at most 

a finite number of vectors / G at which uo ■ k + Vt ■ I \s small. It follows that, having 
fixed an arbitrary constant C, there is at most a finite set of /c's and /'s over which 
\u} ■ k + VL- l\ < C. Denote by S such a set. 
For {k,l) e S consider 

fki{f^) = uj{m) ■ k + Q{m) ■ I 

since f^i is an analytic function it has only isolated zeros. So at most finitely many of 
them fall in [0, L]. The set A is the union over k,l E S of such points. Fix m G [0, i^]\A. 
□ 

Then one can put the system in seminormal form. The explicit computation was 
essentially done in [KP96] (see also [Pos96b,GY03]), obtaining that the matrixes A and 
B are given by 



A 



6 



TT 



U110J2 

3 



4 



B 



6 

TT 



4 



4 



■ (47) 



13 



Remark that, defining the matrixes 5*1 := diag(u;i, ...,ct;„) and 5*2 := diag(a;„+i, a;„+2, ■■■) 
one can write A = ^S^^AS^^ , B = ^S^^BS^^ with 



/3 4 ... 4\ 



4 3 ... 4 



V 4 4 ... 3 / 



B= \ 4 ... 4 



(48) 



With these expressions at hand it is immediate to verify that dctA ^ 0. For what pertains 
Clj by exactiy the same argument in the proof of lemma 4. 1 one has that they are different 
from zero except for at most finitely many values of m e [0, L] . 

Thus, provided m does not belong to a finite subset of [0,L\, theorem 2.1 and its 
corollary 2.1 apply. 

4.2 A nonlinear Schrodinger equation 

Consider the space as in the previous section. Following Poschel [Pos02] we define a 
smoothing operator as follows. Fix a sequence {pj}j>i with the property 

Vj > 1 Pi 7^ , and \pj\ < Cj-"^^^ , d > 1. (49) 

Consider the even, 27r periodic function p{x) :— J2j Pj cos{jx) and define 

T :H'c^ H'c'^^'^ , Vu:= p*u (50) 

where the star denotes convolution (it is defined first extending the function u to an odd 
27r periodic function). 

Remcirk 4.1 It is easy to see that, expanding u in Fourier series 

= ^^isin(ja;), 

the j-th Fourier coefficient ofTu is proportional to PjZj. 

Consider the Hamiltonian system with Hamiltonian function 

H{u,u)^ \u,\'' + P{\ruf) (51) 

with P an analytic function having a zero of order two at the origin. So the equations 
of motion are 

-m^u,, + T[P'[\Tuf)Ttj}j. (52) 

Inserting the Fourier expansion of u the Hamiltonian takes the form (4) with uj = 
and the vector field fulfills (S) with d given by (49). 
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Then one has that (NR) is here violated. So in principle one cannot expect proposi- 
tion 2.1 to hold for this system. However it is clear that one can use a similar procedure 
to obtain a weaker result in which the function G is in resonant normal form, i.e. it 
contains only monomyals Poisson commuting with Hq. Now, an explicit computation, 
identical to that done in [KP96], shows that the so obtained normal form depends on the 
actions only. So one can proceed in the construction and try to look for resonant tori 
and the corresponding periodic orbits. 

The explicit computation shows that also in this case the matrices A and B have 
the structure A — aSiASi , B — 2aS2BSi with the matrixes A and B still given by (48), 
matrixes 5*1 := diag(pi, ...,Pn) and 5*2 := diag(p„+i, p„+2, ....) and a suitable constant a. 
So the determinant of A is still different from zero. The frequencies Clj have now the 
structure 

^i(p) = f - Pi«(p) > Vj > n + 1 
where a is a function of So, except for exceptional choices of {Pj}j>n+i the 

nondegeneracy conditions are fulfilled and the conclusions of theorem 2.1 and its corollary 
apply. 

5 Appendix: Proof of proposition 2.1 

The idea is to proceed as in the proof of the standard Birkhoff normal form theorem, 
i.e. by successive elimination of the nonresonant monomials. As a variant with respect 
to the standard procedure one does not eliminate terms which are at least cubic in the 
variables z. Remark that the estimates involved in the proofs are much more complicated 
than in the finite dimensional case. 

To start with expand P in Taylor series up to order five: P = P3 + P4 + Ps+higher 
order terms. Then we begin by looking for the transformation simplifying P3. So write 

P3 = Pi + Oipf), 

where P^ is composed by the first three terms of the Taylor expansion of P3 in the 
variables z only (so it contains only terms of degree 0,1 and 2 in such variables). We 
use the Lie transform to eliminate from P^ all the nonresonant terms, i.e. we make a 
canonical transformation which is the time 1 flow of an auxiliary Hamiltonian system 
with a Hamiltonian function x of degree 3. By considering the Taylor expansion of 
at zero one has 

Ho^'^Ho + P^ + {x, Ho} + 0{\\z\\') + 0{\\z\f) . (53) 

One wants to determine x so that P3 + {x,Ho} depends on the actions Izj]"^ only. To 
this end we proceed as usual in the theory of Birkhoff normal form. 

Denote by x — {xi,...,Xn) = {zi,...,Zn) the first n variables and take x to be a 
homogeneous polynomial of degree three. Write 

X= E Xnj.jsj.x^^x^'^^'^" (54) 

Ul| + U2| + b-3| + b4|=3 
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with multindexes Ji,j2ij3ij4- For a multi index ji = (j; i, .... j;,J we used the notation 
\ji\ '■= + - • • \ji,n\ and x^^ := x^'^ . . . xit'" , and similarly for a multi index with infinitely 
many components. So, one has 

{X, Ho}= i • Ui - h) + ^ ■ ih - k)) Xh32my^^^^^^^'^^*- 

Write now 

Iii|+li2|+b3|+b4|=3 

and remark that the indexes are here subjected to the further limitation + < 2. 
So one is led to the choice 

■ a. ^'-^^+^-..^0 (56) 

and zero otherwise. Remark that, due to the assumption (NR) the denominators appear- 
ing in the above expression are all different from zero. Moreover, since due to the growth 
of the frequencies they are actually bounded away from zero. Then in order to conclude 
the proof (at least for what concerns the elimination of the third order part) one has to 
ensure that the function x is well defined and that it has a smooth Hamiltonian vector 
field. The terms of x of different degree in z have to be treated in a different way, so 
we will denote by xo, Xi, X2 the homogeneous parts of degree 0,1 and 2 respectively with 
respect to the variables z. 
We need a few lemmas. 

Lemma 5.1 Let R" 3 x ^ f{x) e £^ be a homogeneous bounded polynomial of degree 
r. Write 

jeN»,|j|=r-fe>l 

where is the standard basis of i^. Let {pj^kj'j^^n be a sequence with the property 
\Pjk\ > C and define a function g by 

9{x)= E Y—^'^k. (57) 

jgN",|i|=rfc>l Pi'^ 

Then there exists C such that \\g{x)\\ < CWxW^ . 

Proof. Write g{x) — Y.j QjX^ and remark that the cardinahty of the set over which the 
sum is carried out is finite. We estimate each of the vectors g'j's. One has 

Now the norms of the vectors fj are bounded and therefore the thesis follows. □ 
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Remark 5.1 By the same proof, the same result holds if the space P is substituted by 
the spaces Ti""'^. 

Lemma 5.2 Let R" x 3 {x, z) i— > f{x, z) be a homogeneous bounded polynomial 
of degree r in x, linear and bounded in z. Write 



k>l 
jeN",|j|=r- 

Let {Pj.fel^gNn be as above and define a function g by 

gix,z)^ ^ ^x^Zk. (58) 

k>i Pjk 
3eN",|j|=T- 

Then there exists C such that \g{x,z)\ < C||x||^||2;||. 

Proof. Just write g{x, z) = J2j qA^)^^ ■ Fix 3 ^^'^ study the linear functional gj{z), one 
has 



k>l Pj'^ 



<ll/.lllk/p|| 



where fj is defines in analogy to gj, its norm is the norm as a linear functional, and z/ p 
is the vector of with k-ih. component equal to Zk/ Pjk- Prom this inequality, summing 
over j, the thesis follows. □ 

In order to estimate the vector field of X2 wc will need the following lemma: 

Lemma 5.3 [Lemma A.l of [Pds96a]]. If A — {A^i) is a bounded linear operator on , 
then also B = (Bki) with 

is a bounded linear operator on . 
For the proof we refer to [Pos96a]. 

Lemma 5.4 Let R" x 9 (x, z) ^ f{x, z) & i'^ be a homogeneous bounded polynomial 
of degree r in x linear and bounded in z. Write 

f{x,z)= fjklX^ZkCi. 

k.l>l 
jeN",|j|=r 

Let {Pj,k,i}%^l be a sequence fulfilling 

\Pjki\>Ci{l + \k-l\) (60) 
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and define a function g by 



g{x,z)= E —^'^k. (61) 

k,l>l Pjkl 



Then there exists C such that \\g{x^z)\\ < C||a;||''||2;||. 

Proof. Write g{x,z) = Y,jgj{z)x^ . Fix j and apply lemma 5.3 to such operators ob- 
taining the result. □ 

Remark 5.2 An identical statement holds for functions from R x Tia^s to Tia^s+d- To 
obtain the proof just remark that the boundedness of a linear operator B = {Bki) (gj 
in the proof) as an operator from H""'^ to 7^"''^+'^ is equivalent to the boundedness of 
B :— (vkBkiSi) as an operator from to itself where Vk,si are suitable weights. 

With the above lemmas at hand it easy to estimate the vector field of x- We treat 
exphctly only xi- 

Lemma 5.5 Let Xi be the component linear in z and z of the function x defined by 
(56). Then there exists a constant C such that its vector field is bounded by 

\\X^A^,z)\l,s+d<C\\z\\l,. 

Proof. Write Xi ^ follows 

Xi{x,x, z,z) = {xoi{x,x);z)^2 + (^Xio{x,x);z)^^ . 
Consider the first term. Separating the x, x and z components, its vector field is given 

by 



,-i{^^;z) ,-iXoi[x,x 
dx \ dx I ^2 



Explicitly xm is given by 

\nVHH\=^ l(^ ■ (Jl - J2) + ^^/) 

i>n+l 

It follows that each of the x (and x) components of the vector field has the structure 
considered in lemma 5.2, which therefore gives the estimate of such part of the vector 
field. Concerning the z component lemma 5.1 applies and gives the result. The remaining 
components can be treated exactly in the same way. □ 

The estimate of the vector fields of xo and X2 are obtained in a similar way. In order 
to apply lemma 5.4 to the estimate of the vector field of Xi o^ie has just to remark that 
from (A) and (NR) one has the estimate 



|a;-A; + %-n/| >C(l + |j -/|) . 



Thus we have the following 
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Proposition 5.1 The vector field of the function x defined by (56) fulfills the inequality 

\\X^{z,z)\\a,s+d<C\\z\\l, . 

Then by standard existence and uniqueness theory one has that such vector fields 
defines a unique smooth time 1 flow in a neighbourhood of the origin both in Va,s and in 
Va,s+d- It follows that the transformation is well defined. Transforming the vector field 
of H one gets a vector field having the same smoothness properties of the original one. 
Thus one can iterate the construction and eliminate also the unwanted terms of degree 
four and five. This concludes the proof of proposition 2.1. 
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